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Abstract 

Nonperturbative and perturbative interaction mechanisms of Wilson loops are 
studied within the background field formalism. The first one operates when distance 
between minimal surfaces of the loops is small and may be important for sea quark 
effects and strong decay processes. The second mechanism - perturbative interac- 
tion in nonperturbative confining background is found to be physically dominant for 
all loop configurations characteristic of scattering process. It reduces to perturba- 
tive gluon exchanges at small distances, while at larger distances it corresponds to 
the i-channel exchange of glueball states. Comparison to other approaches is made 
and possible physical applications are discussed. 



1 Introduction 



Interaction of Wilson loops in QCD is the basic element of many physical applications. 
One can mention hadron-hadron scattering amplitude, in particular the phenomenology 
of Pomeron exchange; the assumed color transparency phenomenon; strong hadron decays 
and OZI-forbidden processes etc. There is also considerable interest in calculating Wilson 
loop and Polyakov loop correlators per se, not only in QCD but also in other field theories, 
in particular in supersymmetric Yang-Mills theory. In all cases one starts with connected 
average of two (or more) Wilson loops and tries to calculate it in a kinematic region of 
interest using appropriate field-theoretical technique. It is the aim of the present paper 
to do it in the framework of field correlator method in gluodynamics (see, e.g. review 
and references therein) incorporating both perturbative and nonperturbative contribu- 
tions. Let us briefly remind the basic ideas behind the method. The general Wilson loop 
approach was introduced originally for heavy quarkonia [Q]. For quarks of finite mass one 
can use Feynman-Schwinger representation (see f3[ for a review and references therein) to 
write the meson Green's function as an integral over all possible Wilson loops, formed by 
the quark trajectories and finally to express the meson (and baryon) dynamics in terms 
of gauge-invariant correlators of the field strengths, characterizing the properties of con- 
fining background. When going to the hadron-hadron scattering one can adopt the same 
formalism to express scattering amplitude through the vacuum average of the product 
of two Wilson loops, with subsequent integration over all ensemble of loops. In doing so 
we use the background field formalism to separate nonperturbative gluon configurations 
from perturbative (sometimes called "valence") gluons. Thus the answer will contain two 
parts: purely nonperturbative and perturbative inside noperturbative background, i.e. 
glueball exchanges between Wilson loops. 

For the former it is convenient to use nonabelian Stokes theorem and express the 
answer in terms of the gauge-invariant field correlators and finally via the string tension. 
Perturbative part in the nonperturbative background corresponds to exchange of glueball 
states between loops. We shall keep number of colors N c as a free parameter in what 
follows. It will be argued that the leading term for a typical kinematics of the scattering 
process is the (background-modified) perturbative one. This is in line with the old ob- 
servation that high-energy scattering amplitude is dominated by the Pomeron exchange. 
We do not consider here the leading in 1/N C terms of ordinary Reggeon exchanges, which 
formally refer to one-loop case and are subleading in the high energy limit. We are also 
not discussing pion exchanges which may give the main contribution in some cases at not 
large energies, and will concentrate our attention on the case of theory without dynamical 
quarks, i.e. gluodynamics. In another physical situation, e.g. when accounting for the 
sea-quark loops or for a decay transition of a hadron state, the roles of perturbative and 
nonperturbative mechanisms may change depending on the hadron quantum numbers. 
In all cases however the loop-loop interaction is the starting point of field-correlator for- 
malism application for scattering, strong decay etc. In section 2 we introduce the general 
background formalism for the interaction of Wilson loops. In section 3 the nonperturba- 



2 



tive mechanism is studied in detail. In section 4 perturbative gluon exchange is shown to 
transform into glueball exchange mechanism at large distances. Section 5 is devoted to a 
physical discussion of results and brief comparison with the existing models. 



2 Interaction of Wilson loops in the background field 
formalism 



In this section we are going to exploit the background field formalism 0] in the form 
worked out in |]5[. We refer the interested reader to the cited papers for all the details 
and recall basic steps only briefly. We start with decomposing of gluon field A^(x) into 
nonperturbative background -£> M (x) and perturbative part a M , propagating in the back- 
ground: 

A^ = B^ + (1) 
Total gauge transformation is decomposed as 

B, L -> f/ f [b^ - ~dA U ■ U^a^U (2) 

The principle of separation is of no importance at the level of partition function due to 
obvious identity 

Z = ^J VA„exp(-S[A}) ^^JvB^J Va^ex P (-S[B + a]) (3) 

(here gauge-fixing and ghost terms are assumed to be included into the measure of inte- 
gration). The Wilson loop depends on both and a M : 

W{C) = ± Tr Pexp [ig J (B; + aftfdz^ = 

i M 

= — lirn^Tr J] P (l + i9(B,(z [m] ) + a p (z^))Az^) (4) 

c m=l 

The trace in fundamental representation is normalized as 

Tr i = N c ; Tr t a t b = -6°* 

Our general strategy is the following ||: we expand the correlators under study in 
powers of the field a M , while account for effects caused by the nonperturbative background 
exactly (i.e. without expansion in powers of B^). Namely, one has 

W(C) = W {0) (C) + W {1 \C) + W {2 \C) + ... (5) 
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where, e.g. 

/ \ 

W 0-)(C) = ^ Tr P zu J ait (z)dZp exp ig J B v {u)du v (6) 

c \ c, J 

while W^°'(C) contains only the field and ordering operator P zu takes care of ordering 
a M (z) and B^u). 

Let us define now connected average of two Wilson loops as 

x(Ci,C 2 ) = (^(Ci)iy(C 2 )) - (W(C 1 ))(M/(C 2 )) (7) 
This average can also be expanded in powers of ga M : 

X {Ci,C 2 ) = X (0 \C U C 2 ) + X (»(C 1 ,C 2 ) + X «>(C7 1 ,C 2 ) + ... (8) 

Here X (^1^2) is purely nonperturbative interaction of two Wilson loops and depends 
only on fields B^, while higher terms X^iPii^^) are proportional to the average of 
(ga^) n . One immediately notices that since (W^) is identically zero, the term x (Cx, C 2 ) 
vanishes and the expansion starts with the two-gluon exchange term X ^(Ci, C 2 ). 

In some cases also C-odd exchange contribution (odderon-type) is important, it is 
contained in x^(Ci,C 2 ). In what follows we discuss mostly the purely nonperturbative 
term x' ' (Ci > C2) , and in the last part of the paper also two-gluon exchange x^(C\,C 2 ) 
in section 4. 



3 Nonperturbative interaction of Wilson loops 



We consider in this section the first term in the expansion (H), namely x {Ci,C 2 ) and 
use the contour gauge M to write down the Wilson loop as surface integral 



W{C) = ^TrPexp ^ig J F^(u,2 )*vH) 



(9) 



We have defined in @) F^ u (u, xo) = $ l X0U F IJiU (u)^ l uxq where phase factors along the curve 
L XO u with the edge points x and u are given by 



-• X l\U 



PexpfigJ ? B^dzfA (10) 



Consider Wilson loops defined for two contours C\ and C 2 , where individual minimal 
surfaces will be denoted as S™ m and S™ in , respectively throughout the paper |] The typ- 
ical problem in the discussed framework is to choose optimal integration surfaces in the 

^or any surface we use one and the same letter S for a surface as geometrical object and for its area. 
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integrals of the form (||). The Wilson loop @ is gauge-invariant and surface-independent. 
In the confining regime it obeys the area law, which means that a single surface is unam- 
biguously chosen dynamically (modulo surface fluctuations) and it is naturally assumed to 
be the minimal surface. In the approach proposed in [|7| one can visualize how the minimal 
surface results from the balance of action of all correlators, and higher correlators play 
important role in making the surface smooth and minimal. After this goal is achieved, the 
relative importance of correlators on the minimal surface obeys a strict hierarchy. For a 
single Wilson loop it can be strongly argued (see discussion and references in 0, |§), that 
in case of minimal surface (which is obviously the distinguished surface for given contour) 
the dominant nonperturbative contribution in the cluster expansion of @ comes from 
the lowest, Gaussian correlator of the field strength operators. This property is known as 
Gaussian dominance and it plays important role in all phenomenological applications 
of field correlator method. 

In the physical picture described above appearance of the minimal surface has been 
in some sense a result of field correlator dynamics. However, one can take another view 
by saying that Gaussian dominance corresponds to such profile of the confining string 
worldsheet, which minimizes the total energy of the system. In other words, to calculate 
(W(C)) one may proceed as follows: first, to find the minimal energy of confining string 
configuration (trivially corresponding to minimal area surface in case of single static loop) 
and, as the next step, to calculate average (|) with Gaussian ensemble of correlators 
integrated over this surface. It is this principle which we shall use in what follows to 
choose the shape of the surfaces entering our problem. To illustrate the physics behind 
it, let us consider gauge-invariant Green's function of scalar particle of the mass Mi and 
antiparticle of the mass M 2 interacting with nonabelian gauge field A a ^ (see, e.g. ||) 

oo oo 

= fda 1 J ds 2 J (Vz) xy J {Vz) m exp(-K )(W(C xym )) (11) 



where 

K„ = «?., + M is2 + I J* J d r 2 (g 2 (12) 

and the contour C xyyx is formed by dynamical trajectories z, z and the lines L xx , L yy . 
Green's function (|TTD encodes all information about the spectrum of the system. Of 
special interest are particular cases when some trajectories in the integrals in (|ll|) are 
singled out kinematically as providing the main contribution. Phenomenologically the 
most important cases are represented by the large mass (nonrelativistic) limit, where 
the particle trajectories are close to parallel to the temporal (Euclidean) axis, and large- 
momentum limit, where eikonal approximation is applicable. Let us consider the former 
case in more details and assume for simplicity that X4 = £4, = f/4, y± — X4 = T, and 
center-of-mass motion of the considered two-body system is absent. The Green's function 
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(PI) can be spectrally decomposed as 

G® (x, x; y,y) = Yl ^{n}^ ~ ^{n}(y ~ v) exp{-E™T) (13) 

{n} 

where T stays for the Euclidean time interval while ip{ n }, denote wave function and 
energy of the state with quantum numbers {n}, respectively. Comparing ([□]) and flT3| ) 
we notice that the problem of minimizing energy levels is in close correspondence with 
the problem of maximizing average value of the Wilson loop (W{C)) for given contour, 
prescribed by kinematics. This correspondence becomes straightforward in classical limit 
of static rectangular loop where area law dictates 

(W(R x T)) » exp(-ET) 

and the static energy is given by E = aR, which corresponds to straight line string and 
minimal area. Notice that we confine ourselves to the leading area law terms in this paper 
and do not consider perimeter contributions. 

For two or more Wilson loops the situation is more complicated, but general frame- 
work is the same. Indeed, one can consider multi-point Green's functions, e.g. two-meson 
one 

g [4] = (*L*L*^**-> ( 14 ) 

where 

and again adopt Feynman-Schwinger representation to get 

G [4] ~ [[***]- {W(C 1 )W(C 2 )) + ... (15) 

where stars [* * *] stay for integration measure while dots denote other possible ways to 
contract field operators in the definition ( |14D and get, correspondingly, other geometries 
of contours C\, C^. On the other hand, the Green's function ( |T4] ) describing interacting 
two-meson system, must be dominated at large separation/time by the confining strings 
configuration of lowest total energy, in direct analogy with one meson case. 

It is physically obvious that for two well separated loops the minimal configura- 
tion of confining strings is given by two individual minimal surfaces with no common 
points of intersection. From the point of view of the expression (pl|) it means that aver- 
age (W(Ci)W(C2)) factorizes into (W(Ci)) (PF(C 2 )} and no nonperturbative interaction 
takes place. Notice that perturbative interaction is mediated in our picture by propa- 
gating degrees of freedom, while could receive nonzero contributions only because 
of rearrangement of confining strings with respect to the free case, if it happens. In the 
approach developed in the present paper it is always assumed that the gluon correlation 
length T g (which also defines the radius of the string) is taken vanishingly small, while 
the string tension a is kept fixed and therefore is expressed in terms of a alone. This 
is in contrast to another approach |J, to be discussed later in the paper. 
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Consider now a process of decreasing the distance between the loops. For classical 
films, as it is easily checked in experiments with soap films, at some point there appears 
the situation, when the minimal surface, bounded by contours C\ and C2, which we shall 
call Su, becomes smaller, than the sum of S™ in and S™ tn . At this point, according to 
minimum energy principle, there should appear a common surface Su which constitutes 
the interaction between the loops. Indeed the classical films (e.g. soap films) satisfy the 
stated above principle. For the simplest case of two concentric circles of equal radius R 
at distance h it was Leonard Euler who demonstrated that the corresponding minimal 
surface Su is a catenoid and it exists for any h such that the equation 

cosh ( -J^—] (16) 



Rmin \2-Rj7 

has a solution. Here R m in defines the so called minimal radius of the catenoid. At some 
critical distance h cr the solution ceases to exist and the minimal surface coincides with 
two disconnected disks.Q 

The case of QCD is more complicated, first of all because of nonabelian nature of 
the gluon field. In particular, it is possible strictly speaking to take in ([]) integration 
surface with the topology different from that of the disc only at the price of introducing 
additional holonomy factors.^ It is also worth mentioning that the surface corresponding 
to the global minimum of energy can be glued with itself in nontrivial way, as it happens, 
for example, in baryon correlators. In fact, soap films demonstrate a similar behavior 
to some extent: they easily form "string junction" with three worldsheets having one 
common intersection line and forming e.g. "Mercedes-star" configurations. Moreover, 
contrary to the case of one Wilson loop there is no, in some sense, a unique distinguished 
string worldsheet configuration for the product of two or more loops. Let us illustrate 
this on the example of two coinciding and oppositely oriented Wilson loops, where the 
following well known identity holds true 

\W fund \ 2 = ± + (1 - ^ • W adj (17) 

where Wf un d is given just by (|), while W a dj corresponds to (f|) after replacement of fun- 
damental SU(N C ) generators t a with generators [T a ]b c = —if abc in adjoint representation 
and proper change of trace normalization.^ The terms in the r.h.s. of (|17D are in one- 
to-one correspondence with decomposition of the product of two fundamental SU(N C ) 
representations into singlet and adjoint: N c ® N c = 1 © — 1). Accordingly, the answer 

2 Actually the analog of condition ( pi^ ) is violated earlier, at some distance h cr < h cr . The catenoid 
solution is locally stable but globally unstable for h cr < h < h cr . At h > h cr local stability is lost. We 
do not make a distinction between h cr and h cr in our consideration and speaking about critical distance 
we always have in mind minimal energy condition. 

3 See, e.g. It is in contrast with abelian Stokes theorem where one can easily integrate two-form 
field strength over noncontractible surface. 

4 Traces (both fundamental and adjoint) are normalized to unity in (fL7ft. 
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is given by superposition of the state where two fundamental confining strings are totally 
annihilated (corresponding to the unit representation and constant term in (0)), and the 
state where confining strings are summed into one adjoint string (described by the second 
term in the r.h.s. of (|I7D). It is clear that for large contours the second contribution dies 
out and the minimal energy state is given by the singlet state. 

Till the present moment we have discussed minimization of energies, paying no at- 
tention to pre-exponential factors like 1/N% in (|T7[) . The latter can become important in 
studies of large- N c behavior of the correlators. We shall come back to this question later 
on. 

The physical picture outlined above is in fact quite general. Therefore we adopt 
the following algorithm for computing nonperturbative Wilson loop correlator. First, 
for given geometry of the contours, we have to find confining string configuration which 
gives dominant contribution to the quantity (W(Ci)W(C 2 )) (notice, that there might 
be different surfaces for different representations, which appear in the expansion of the 
product of fundamental Wilson loops, as in (17) ). If the configuration, providing the 
maximum coincides with individual minimal surfaces, we conclude that nonperturbative 
interaction is absent. As the second step, we calculate correlator x ^ {Ci,C 2 ) via (||) 
in terms of gauge-invariant field correlators. We assume that the ensemble of correlators 
obeys the same hierarchy on this surface as it does on the minimal surface for single Wilson 
loop, i.e. exhibits Gaussian dominance. The strategy adopted below for computation of 
higher terms in the expansion (|8D is slightly different but basic physical arguments for 
choosing minimal surfaces remain the same (see section 4). 

We are keeping generality at the moment and perform the vacuum averaging of 
Wilson loops, i.e. we explicitly calculate x {Ci,C 2 ). Suppose that surfaces Si , S 2 have 
already been chosen according to our criteria. For field-strength tensors, belonging to 
surfaces Si and S 2 and gauge-transported to the same point x we define: 



(Fda 



■it 



(u,x )da liu (u) ; u G S x 



(Fda) ( \v) = F IJ , I/ (v,Xo)da lxu (v) ; v e S 2 
One can now write down the product of two Wilson loops in matrix form as 

N 2 C {W{C X )W{C 2 )) = 



;is) 



Tr iTr 2 P 



12 









( 


exp 


^ ig f(Fdaf\u)\ 




exp 








\ 



ig i(Fda) {2 \v) 



(19) 



012P2 



where the traces Tr 1 and Tr 2 go over indices carrying subscripts 1 and 2 respectively and 
the ordering operator P12 orders the products of matrices in a proper way according to 
the definition of P-exponent (0). 
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One can derive the following rule for vacuum averaging of several matrix operators, 
transported to one point (it is easy to show that these relations are gauge-invariant under 
field-independent gauge transformations) 

([F(uM, x o )..F(u^,x o )] ai01 [F(v^\x )..F(v {m \ x )} a ^ 2 ) = 



^ai/3 2 ^a 2 /3i 

m- 1 



(Tr (F(uW, x )..F(u^,x )F{v^,x )..F(v^\ x ))) 



-— (Tr (F(uU, x )..F(u (k \ x ))Tr x )..F(v (m \ x ))) 



(Tr(F(u^\x )..F(u^\x ))Tr{F{v^,x )..F(v^,x ))) 



-—(Tr (F(uW,x Q )..F(u( k \ x Q )F(v^, x )..F(v^\x ))} 
For lowest Gaussian correlator when k = m = 1 one has from (2D) 



{[F(U, So)]ai/9i [F(V, X )} a2 /3 2 ) 



(Tr F(u,x )F(v,x )) 
N?-l 



Nr. 



aif3i°a 2 (3 2 



(20) 



(21) 



The above relations are valid in theories without global color symmetry breaking. Now 
one can proceed with matrix cumulant expansion (see, e.g. [ljj): 



N 2 C {W{C X )W{C 2 )) = Tr xTr 2 P 12 exp £ 



{ig) n 



\n=l 



ni 



«F(l)...F(n)» 



(221 



where we use the notation 



F(k) = I (Fda) w (u^ k) ) + I (Fda) {2) (v w ) 



(23) 



Si 



s 2 



Here double brackets ((...)) denote irreducible correlators (see definition in [|I|, |TTJ). The 
quantity F(k) carries four independent fundamenatal color indices and the traces Tr 1; Tr 2 
go over indices corresponding to F^ and F^ 2 \ respectively. 



Expression (J22|) provides the basis for our discussion. For a single loop the corre- 
sponding cluster expansion is given by 



(W(C 1 )) = — TnPexp 



E ((/ •••••/ (Fda)^(u^))) I (24) 
Si Si / 



and in area law regime (which means that typical sizes of the loops are larger than the 
gluon correlation length T g ) one gets 



{W(d)) oc exp(-aS™ n ) 



(25) 
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The string tension a is given by 

aSr n = d <u {1) ) J da( U W)(Tr gF(u^ ,x Q )gF(u^ ,x )) + ... (26) 

C 5i Si 

where the dots denote higher non-Gaussian terms. As it was mentioned above, we do not 
take into account perimeter terms, it is implicitly supposed that all considered loops are 
large enough in this sense. It is straightforward to rewrite ( [22] ) combined with (|23p and 
( p0|) in the following way 

Nl(W{C x )W{C*)) = Tr iTr 2 exp (l • (A + Ai) + e • A e ) (27) 

where 

A = — (tS\ - aS 2 (28) 

and terms Ai and A e contain correlators of powers of F defined on different surfaces and 
hence provide contribution to ■ They are as follows: 

Ai = ± N( S_ 1) J M^) J da(v^){TrF(u^,x )F(v^,x )) + ... (29) 
c ^ c 'si s 2 

and 

2 

Ae = T WZT[ J da ^ J da(v {1) )(Tr F(u^,x )F(v^,x )) + ... (30) 

c Sx S 2 

where the dots again denote higher non-Gaussian terms. The upper (lower) sign in the 
above expressions corresponds to the case of parallel (opposite) orientation of the contours 
G\ and C 2 and surfaces S\ and S 2 (we assume that the orientation of the surface is fixed 
by the orientation of its boundary contour). 

The matrix structures 1 and e introduced above are given in index notation by 

We are to examine what algebra the matrices 1 and e do obey. In the chosen approxima- 
tion it encodes the effect of ordering, performed by the operator P12. Consider the case 
of parallel orientation. In this case ordering of the fields F^ 1 ' 2 ^ in the Wilson loop product 
coincide for both loops, which for the matrices 1 and e mean, that they should be multi- 
plied from the right with respect to indices carrying subscript " 1" and corresponding to 
the surface Si, and also with respect to indices carrying subscript "2" and corresponding 
to the surface S2. It is easy to check, that this requirement transforms into the following 
relations: 

e • i = i • e = e ; P = i ; e 2 = l (32) 

In case of antiparallel orientations of surfaces the matrix e should be multiplied from 
the right with respect to indices carrying subscript "2" but from the left with respect to 
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indices carrying subscript " 1" , which corresponds to the fact of opposite ordering of the 
fields and F^. It results in the algebra different from (|32|): 



e ■ 1 = 1 • e = e ; 1=1; e = N c e 
We also notice that in both cases Tr iTr 2 1 = iV? and Tr iTr 2 e = AT C . 



(33) 



With (0), fl33|) at hand, we can easily compute (|27|) . For parallel orientation of 
surfaces one gets 



x (0) (Ci>c 2 ) 



K 1_ iv) 6XP ( A + A l" A e ) + 



+ 2 (1 



A^ r 



exp (A + Ax + A e 



- exp (A™ n ) 



(34) 



where the last term corresponds to the product of averages of two loops. For oppositely 
directed contours the result is 



exp (A + Ai + N C A C 



1-^Jexp (Ao + Ai) 



exp (AZ n ' 



J 

(35) 



Expressions (|34|) , (|35|) together with the prescription for the choice of optimal inte- 
gration surface provide the answer for nonperturbative interaction term. Let us come to 
concrete examples and first consider the simplest possible case when Si = S™ m , S 2 = S 2 nm 
and S! 2 nm C S™ m (see Fig.l). As is was explained in details above, we always assume 
Gaussian dominance. One easily gets the following result for the contours with linear 
sizes greater than T g and omitting perimeter contributions: 



Ai 



and @§, (|3|) become 



"A? - 1 



oS™ m 



X (0) (Ci, C 2 ) = exp(-o-S , r n - aSt n 

2a S^ 



A -i- ^ c nmin 



1 / 1\ /2cx5™ n \ 

2 {'-nJ^ [n^T + 



(36) 



K i+ ^ ,<,xi> 



N c + 1 



— exp( — & Si — ab. 



min\ 
2 J 



(37) 



for parallel orientations and 



+ 1-^ exp 



2aS^ in 



-lexp(2a^ n ) + 

. c 

~ omin _ omin\ 

exp{—ao 1 — cro 2 )■ 



(38) 
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for opposite orientations. It is worth reminding that S™ m > S™ m by our convention. 

Notice that in case of coinciding, but oppositely directed contours (i.e. C — C\ — 
[C2P), expression (|3^) reproduces (0), if adjoint string tension is given by the Casimir 
ratio 

2N 2 

Vadj = JJT^T a ( 39 ) 

c 

This result is to be expected since, as it was already mentioned, Gaussian dominance 
yields Casimir scaling |J. If both contours C\, C2 lie on the same plane and C2 is inside 
Ci, the geometry becomes effectively two dimensional and the results (|37|) , (|38| ) coincide 



with formulas obtained in slightly different way in [12]. The same expressions hold true in 
1+1 dimensional Yang-Mills theory, where one has just two dimensional geometry and, on 
the other hand, exact Gaussian picture. One also look upon (0), ( p8[ ) as an algebraic rule 
of adding up parallel or antiparallel fundamental fluxes which illustrate the decomposition 
3® 3 = 3 ©6 and 3® 3 = 1©8 respectively, with the string tension in each representation 
given by the Casimir scaling law . 

We are now interested in the case of contours separated by distances greater than 
T g . By way of example let us calculate purely nonperturbative correlator of two Wilson 
loops {W {C\)W (C2)) for simple rectangular geometry of the contours. We choose two 
rectangular contours RxT lying on parallel planes, at distance h from each other (see Figs. 
2, 3). We suppose that T 3> R and will not take care of subleading 1/T terms. If h is of 
the order of T g , one comes back to the case described by (0), (|38|). We take the distance 
h such that R > h > T g , where the nonperturbative regime is supposed to play a role. 
One could still choose for Si and 5*2 the corresponding minimal surfaces S™ m and S 1 ™". 
In this case the correlator would be equal to zero up to exponentially small terms of 
the kind exp(—2h/T g ), as it is clearly seen from (p9|), (|30|). Such configuration of surfaces 
however does not correspond to the minimal energy condition formulated above and is 
therefore unstable. It is clear that the subtraction term, the last exponent on the r.h.s. of 
(j37[), ( p8|) stays intact, since it is a product of single loop averages. Other terms, on the 
contrary, strongly depend on the profile of the strings, which defines the integration surface 
and is to be chosen according to dynamical minimal energy condition. Consider first the 
case of opposite loop orientations and let us examine different choices of surfaces. As a 
trivial example we might adopt the same choice as above, namely Si = S™ m ; S2 = S™ m . 
The first term in square brackets in (^) contributes to (W \C\)W (C2)) as 

-Uxp(-2ai?T) 

c 

as it should be since A = A™ in = aS?™ + aS^ in = loRT and Ai = A e = for this 
choice of surfaces (up to exponentially small terms ~ exp(—2h/T g ), which we always omit 
in the paper). Correspondingly, the second term in square brackets in (|35| ) would be 

(l-±)e^(-2aRT) 
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with the sum of two giving expected answer (W (Ci)) (W (C2)) = exp(— 2aRT) and hence 
= 0. One immediately sees that such choice is not optimal if h is small. Instead, if 
we choose Si as a minimal enveloping surface with boundary on Ci and coinciding with 
S 2 = S™ tn inside C 2 (due to apparent symmetry of our problem, one could of course easily 
interchange indices '1' and '2'), we get for the first term in square brackets in (|35| ) 

— exp (-2ahT) 

c 

where 2hT = Si 2 — Si — S 2 . The second term contributes to (W(Ci)W(C 2 )) as 

^1 - -^j exp (-2ahT - a adj RT) 

where a a dj is given by ( j3l?D and this contribution is always subleading with respect to the 
former one if R 3> h and N c is not exponentially large (see below). So these two different 
choices give different answers for Wilson loop correlator: 

(W(Ci)W(C 2 )} = exp(-2aRT) (40) 

in the first case (two individual minimal surfaces) and 

(W(Ci)W(C 2 )) = i exp (-2ahT) + (l - ±j exp {-2ahT - a adj RT) (41) 

in the second case (enveloping geometry). According to our criteria, the answer which is 
dominant should be chosen as correct, since it corresponds to actual string configuration. 
It is seen, that there is a critical distance between loops 

h crit ^R-^-\ogN c (42) 
oT 

in our problem.^ For h < h cr u confining strings rearrange themselves^ with respect 
to noninteracting case, which is encoded in expression (|4T|). Correspondingly, one has 
nonzero x^- F° r larger h they do not interact and x {Ci,C 2 ) vanishes. It is important 
that in our picture it happens dynamically, in particular, one cannot just naively take 
large iV c -limit in (|34D , (|35D . On the other hand, if h is kept fixed, then it is clear from 
([42]) that in large iV c limit x^ should vanish. It is worth saying a few words about the 
meaning of T in (flUP, (f4l"D . It is seen that with T going to infinity, h crit is increasing 



and approaching R. Therefore for long static loops minimal energy state will always 
eventually win, as it should be. If however one studies not correlators of static loops 
but some physical process like hadron scattering, pre-exponent factors are important and 
T corresponds to the typical interaction time. In particular cases it can be quite small. 



5 For abelian confining strings including soap films the last term in (|4^) is absent. 

6 See also [O, where static multiquark interactions were studied in strong coupling expansion regime. 
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Expression (|4^) then physically means that if T is small, strings does not have enough 
time to rearrange and nonperturbative interaction does not take place. 



To summarize, the answer is given by 

X (0) = ^ exp (-2ahT) + (l - ^ exp (-2ahT - a adj RT) - exp {-2a RT) (43) 

if h < h cr it and x^ = if /i > h CT n. We now come to the case of parallel orientations. 
Let us consider two geometries analyzed above. The case of individual minimal surfaces 
is identical for parallel and opposite orientations and has just been considered. The 
enveloping geometry gives for the first term in square brackets in 

while the second term becomes 

2aRT(N c + 2)\ 



K l + w) exp (- 2 " hT 



N c + 1 



and is always subleading. In particular case N c = 3 the former term corresponds to the 
creation of the string in 3 representation of the length R which has the same tension 
as the fundamental string while the second term describes sextet representation string 
formation 6 = Symm{3 ® 3} with tension 2(N C + 2)/(N c + 1) = 2.5 times larger than 
fundamental one. In this particular case of N c = 3 one can additionally consider double 
Y-shape profile, shown on Fig. 3, with the result for the first term in square brackets in 

(0) x 

- exp(—V3ahT — aRT) 

Since < 2, this term is seen to be dominant over enveloping geometry in iV c = 3 case. 

As in opposite orientations case there appears critical length h cr u which for N c = 3 
is given by 



The expression for becomes (iV c = 3): 

1 



x (0)(C 1 ,C 2 )=exp(-a 1 S 12 ) 



2 / 5aS^ 



+ - exp 



3 r V 2 



exp (-<rS 3 ) 



exp(-2aRT) (45) 



where S12 represents the boundary surface Si 2 = 4/iT '/ ' \/3 and S3 = (R — h/\/3)T is the 
common part of Si and S 2 . Notice that l cr i t = R — h crit /^/3 > 0, therefore one is never 
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in the situation of S3 shrinking to zero. In some sense l cr u in ( [45]) plays the role of R m %n 
from (Jl6|). 

The physical picture is the same as for opposite orientations - at small distances h, 
hcrit > h > T g preferable string configuration is given by double Y-shape profile shown 
on Fig. 3 while at larger distances there is no common string state formation and hence 
nonperturbative interaction is absent, x (Ci,^) — 0. 

It is important to compare at this point our approach to another approach, which is 
also based on field correlator method, but follows another logic. This well-known strategy 
started in M has many successful phenomenological applications and modifications (see, 



e.g. [0). There are two essential differences to be mentioned. The first one comes 
from the fact that we are working in Euclidean metric rather than in Minkowski. It is 
important to stress, that in our approach the vacuum averaging over gluon fields (which 
eventually creates the string) is to be done in Euclidean space with Euclidean gluon 
configurations, and transition to Minkowski space can be done only after this vacuum 
averaging. This is intimately connected to the realization that vacuum configurations are 
of tunnelling type and it is not legitimate to continue them analytically into Minkowski 
space. However, as soon as the fields are integrated out the transition to the latter is not a 
problem and as we demonstrated in ]n| the QCD string is present also on the light cone. 



Correspondingly, we do not have a notion of light-like surfaces with vanishing areas, which 
is of importance in the formalism developed in ||. The second point is a different physical 
picture behind the purely nonperturbative interaction in our model. The magnitude of 
X^°\Ci, C2) depends in our approach on the chosen string worldsheet configuration which, 
in its turn, is determined by the dynamical minimal energy requirement, discussed above. 
However, we do not study effects of string overlap, caused by finite thickness of the 
strings. If the minimal surface requirement dictates Si^ = S™", which always happen 
for well separated loops, there is no purely nonperturbative interaction in our picture, i.e. 
X^(Ci,C 2 ) = 0. This is in sharp contrast to the approach || where all physical effect 
is due to nonzero value of T g which bring about nonzero overlap between hadron Wilson 
loops. In other words, the expressions ([37]), (38) have nontrivial limit if T g — > 0, while 



string tension a ~ (F )Tg is kept fixed, contrary to the corresponding expressions from 



|9|, [14] which vanish in this limit. Accordingly the hadron-hadron cross-section in that 
approach is proportional to T g in the tenth power |TJJ and is very sensitive to exact value 



of T g . As it is, both approaches are viable in their regions of values of T g , and it is very 
important to extract this value from lattice data, analytically or phenomenologically. The 
existing calculations yield values in the range 0.35 Fm > T g > 0.1 Fm for different lattice 
and analytic procedures (see [|nj and references therein). We assume in our approach that 
T g is closer to lower limit and therefore the string overlap effect proportional to higher 
powers of T g is a small correction to our basic mechanism. One should stress however, 
that if even T g would be large to justify approach of |9|, [L4J] at moderate energies, the 
gluon-exchange-generated Pomeron mechanism is always dominating asymptotically at 
large energies. 
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To conclude this section let us say a few words about the case of distant loops. The 
general expression (|22[) is applicable here as well. In Gaussian approximation one still 
has fl27p. Important change however is that now correlator (F^F^) does not get any 
contributions since the surfaces S\, S2 should be chosen according to our principle, as 
independent and not intersecting minimal surfaces. Hence nonperturbative contribution 
vanishes in this case, up to the string overlap effects mentioned above. 

In other words when contours are distant from each other, it is always preferable to 
deform Su to two discs, corresponding to minimal surfaces for each contour plus whatever 
thin tube connecting these surfaces through the point xq, thus reducing x (CiiCz) to 
zero. This is a sign, that purely nonperturbative contribution vanishes, and one should 
consider next terms in the expansion (JD, namely x (CijC^)- It is important to realize 
that the surface entering ([5]) is not dynamicalQ, which reflects itself in the possibility of 
infinite squeezing of the tube, connecting two distant minimal surfaces. The situation 
changes however when one includes perturbative gluons, propagating inside its wall and 
forming a physical glueball state in this way. It is actually this glueball exchange mech- 
anism, which corresponds to the term \ (Ci)C2) when background field is taken into 
account. We study this term in the next secton. 



4 Glueball exchange interaction 

We now turn to the term x^ (Ci> C2) in the expansion (|8|). We shall see that in the absence 
of NP background this term reduces to the purely perturbative two-gluon exchange term 



suggested in |T7| as a basic element of Pomeron exchange. To be more precise, let us 
consider first the case of no background fields. In the Feynman gauge for a M , one has in 
the lowest order 

K^y)) = ^ a \^ { l_ yY ( 46 ) 

Our expression for x will have the same form as (|27|) , where one should keep only 
nonperturbative background field in diagonal correlators (Tr F 1 F 1 ) while terms pro- 
portional to (TtF^-F 2 ) contain only perturbative exchange. These perturbative exchanges 
are modified however by the presence of nonperturbative background. To take it into ac- 
count, one has to perform averaging in two steps: first in valence (perturbative) field a M 
and second - in background field B^. 



(2) 

V-2 



(W{C X )W{C 2 )) = |^Tr jTr 2 P 12 J dx% J dxg / dyg> / dy 



7 It has come from the Stokes theorem and therefore could be arbitrary, subject to our principle of 
minimal action. 
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• ((a(x^)a(y^))(a(x^)a(y^)} + (a(x^)a(y^))(a(x^)a(y^)}) 
where short-hand notation was used a(x^) = af}(x^ 



(47) 

The coordinates x^ l \ x^ are 
ordered along the contour C\ as well as y^ l \ y^ 2 ' are ordered along the contour Ci- Notice 



also gauge- invariance of fllTD due to transformation law (^j). 

Before proceeding further one is to define the dependence of (a^{x)a v {%j)) on back- 
ground fields. To this end it is convenient to use Feynman-Schwinger representation for 
gluon Green's function || and represent it as 



G% = K(x)al(y)) = J ds(Vz) xy exp(-K )^ u (x,y) 



(4, 



where 



and 



®t(x } y) 



Kn 




f^a exp 



y 



rx V [2g I drF(z{r)) 



and a, b are adjoint color indices, whereas [i, v - Lorentz indices, i.e. [F] 



ab 



ab 



ft]' 



abc 



To understand better the topology of the resulting construction, it is useful to consider 
large N c limit. One can write in this limit for adjoint phase factors in ( (48f) 



[t a } a ^ ab (x iy )[t% s 



a exp 




1 



(49) 



where <& a g(x,y) is parallel transporter in fundamental representation. Expression (]49f) 
exemplifies well known 't Hooft's rule for replacing gluon line by double adjoint in large- 
one obviously obtains two new Wilson loops Cn and 



N c limit. Inserting (f49|) into 
C\2 instead of previous C\ and C 2 : each initial loop is now divided by two gluon emis- 
sions/absorptions into two arcs which are connected by double lines of gluon propagators 
(see Figs. 4,5). At small N c this construction goes over into that of two fundamental 
loops C\ and C 2 connected by two adjoint lines and final result will amount to replacing 
double fundamental string worldsheet by one adjoint string wordsheet. In terms of string 
tensions it correspond to replacement of 2a by 9a/ 4 in Gaussian approximation. We will 
keep large N c limit and replacement fl49|) in what follows. 



The averaging over background fields leads to the following result: 



{{W{C l )W{C 2 ))^) B = X ^\C 1 ,C 2 
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= jp I ds i VzW I ds * Vz{2) expi-K^-K^) (W(C 12 )W(C 12 ')) (50) 

C L 13 L12' 

where g 2 = g 2 N c . Here contours Ci 2 and Gyi comprise pieces of C\ and C 2 connected 
by two double fundamental lines L\i and L12 . It is understood that surfaces S12, 5*12' 
are subjects of our general assumption about minimal action. This gives different forms 
depending on the distance between original loops C\ and C 2 , see below. It is also under- 
stood that gluon spin operators 2gF(z) are to be placed on the gluon trajectories L12 and 
L12' in accordance with (fR|). It will produce gluon spin interaction terms which influence 
glueball Green's function; to simplify discussion we omit these terms at the moment. 

Now we can use large-iV c factorization property for the product (W \C\2)W (Cu ')) 
and use area law asymptotics for each piece, i.e. for surfaces S12 and 5*12'. One obtains 

X (i) (C l ,C 2 ) = -^ J dsxDz® J ds 2 Vz® exp(-^ 1} - K^) exp(-a(S 12 + S l2 ')) (51) 

c £12 W 



To define the profiles of the surfaces we shall use the same principle outlined above, i.e. 
we require an effective value of area (5*12 + Su') averaged over possible gluon trajectories 
L12 and L12 to be minimal. The result will of course strongly depend on relative positions 
and orientations of the contours C\ and C 2 . In the first case when both loops lie on the 
same plane and C 2 is entirely inside C±, it is clear that the sum Su + Su does not depend 
on trajectories L 12 and L 12 ', and one has S12 + S12 = S\ — S%. Thus one obtains effectively 
the surface Si with the hole due to C 2 , i.e. a construction which already has appeared 
in purely nonperturbative term for opposite oriented contours Ci,C 2 , but now with two 
valence gluons connecting contours C\ and C 2 , see Fig. 4. 

However for large enough distances the true minimum of S12 + S12' is reached by 
another construction - when two contours C\ and C 2 are connected by narrow strip 
formed by trajectories L 12 and L 12 ' with the double (adjoint) string worldsheet between 
them (Fig. 5). This narrow strip is nothing but the glueball Green's function and the 
width of the strip is equal to the average size of the lowest mass glueball, i.e. around 
0.5 Fm. Notice that due to kinetic terms in ([51]) this strip is dynamical (contrary to the 
nonperturbative case) and cannot be shrinked. This is a typical construction for high- 
energy scattering amplitude when the glueball exchange diagram is gradually replaced by 
glueball Regge trajectory exchange, i.e. by Pomeron exchange, which persists to larger 
experimentally accessible energies. 

To demonstrate that explicitly, one should rewrite expression ([51]) directly in terms 
of glueball Green's function: 

(W(C 1 )W(C 2 )) 

P 12 / <*$ / dxg / dy% J dy% [Gft£ (xW,xW\yW,yW) + (yW <- y^)] (52) 

Ci Ci C2 C2 
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where G V a^{ x i , y^) is two-gluon glueball Green's function, describing propaga- 

tion from points x^\ to y^ l \y^ 2 \ which has the Feynman-Schwinger representation as 
in (|5"ID . The spectrum of this Green's function (with spin terms included) was computed 



analytically in |18 



In case all points are close to each other, i.e. 

\ X H) _ X U)\ < T g ; \x {i) - y ij) \ < T g 
one can replace glueball Green's function by a product of free gluon propagators 



rpw ( T {1) T (2)| ? ,(l) „(2)\ ^ ( A OfMivAMM ( ro\ 

^MiwA* ' X W >V ) 9 ( 27 r)4(x(l) -y(l))2( X (2) -j/(2))2 ^ 



Another asymptotics is available when both \x^' — y^'\ are large, then the spectral de- 
composition is possible: 



E {v {1 \ y {2) ) ■ exp f-M n 

n \ 



+ x (2) 3,(1) + y (2) 



(54) 



Since the lowest glueball is rather heavy, M ~ 1.5 GeV, one expect fast decrease of 
when distance between loops is growing: 

X {A \h)~exp(-M \h\) 

The situation is qualitatively similar to the one studied in |TP|, pp| , where gluon was 
assumed to have effective mass m g ~ 0.9 GeV. One expects dipole-dipole cross section 
around a few mb in this case, when a s is of the order of one. To obtain realistic large 
hadron-hadron scattering one needs glueball exchange to be reggeized, in which case 
radius of interaction grows logarithmically . In particular case of BFKL Pomeron this 



picture was studied in ||22|| . Our picture differs from that of BFKL, since nonperturbative 
background is taken into account. Let us consider as an example the problem of high- 
energy forward onium-onium scattering. Since the interaction time between particles at 
high energies is much smaller than the typical interaction time for quarks inside onium, 
one can consider the onium in this process as free quark-antiquark pair (see, e.g. |BJ). 
The small radius of the onium compared with typical transversal length scales of the 
problem dictates x dominance over i n the problem since possible nonperturbative 
string configurations which could contribute to x^ are strongly suppressed over individual 
minimal noninteracting strings. 

Since it is more convenient to study the scattering of systems in given quantum 
states rather than scattering of Wilson loops, we are to switch from ( [T5D to spectral 
decomposition of the form (|13D and take only one term, corresponding to the scattering 
of particular states. The resulting expression coincides (up to normalization factor) with 
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the scattering amplitude (see, e.g. f23fl). Since in the Feynman gauge for the field we 
have 

G^ 2 {x^\x^\yV) = 8» lV1 5^\^\xV\yV,yW) (55) 



the answer can be straightforwardly obtained using the same strategy as in |]T9 

1 1 



Tjorw =i J d 2 p x J d 2 p 2 J dzi J dz 2 \ip(pi, z 1 )\ 2 \ip(p 2 , z 2 )\ 2 F(p u p 2 , Q = 0) (56) 

where 



o o 



N — 1 r 

F( Pl ,p 2 ,Q = 0) = ^-^- 2 j d 2 kG^ 2 \k,Q = 0y 

■ (2 — exp(— ikpi) — exp(ifcpi))(2 — exp(— ikp 2 ) — exp(ikp 2 )) (57) 

In the above expression G^(k,Q) is Fourier transform of ( pSS] ) with respect to total mo- 
mentum pi + p 2 = Q and relative momentum pi — p 2 = Q — 2k, the former is equal 
to zero for forward scattering amplitude (we also suppose in (|57|) vanishing transverse 
momenta of onia). The mixed representation wave function ip(pi, zi) defined on the light- 
cone describes the state of color dipole with transverse size pi and fraction of total onium 
light-cone momentum pf carried by quark Z\. We omit spinor indices, assuming proper 
summation over them. If one will "turn off" confinement (i.e. in our formalism put the 
confining background field to zero everywhere), the Green's function in the leading order 
of perturbation theory would be just a product of two gluon propagators 

G fr^Q) = J^ m (58) 

and inserting (^) into (|57|) one returns to the results of [19|. One would expect that 



the effects of confinement suppress amplitude [|19j in two different ways: first, because of 
the mass gap (and actually quite large mass even of the lightest glueball), and, second, 
due to the glueball wave function fast decrease at large relative distances (see (0)). 
These properties solve the artefact of color Van der Waals forces, appearing in purely 
perturbative dipole-dipole interaction. 



5 Conclusions 



In the present paper we discussed interactions of Wilson loops in confining theory, hav- 
ing in mind gluodynamics as a concrete example. The effects of confinement were taken 
into account in the formalism of perturbation theory in confining background. We have 
described the background by gauge-invariant Gaussian correlator with small correlation 
length, which is supported by lattice and analytic calculations. Two main physically dif- 
ferent mechanisms of interaction were analyzed. The first one, which we call nonpertur- 
bative, refers to the process of confining string rearrangement, which can be energetically 
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preferable for particular geometries of the contours. In this way a common surface of two 
contours Ci, C 2 is created and in case of opposite orientation this surface is a ring between 
C\ and C*2 (with a hole inside the smaller loop). This mechanism has a direct classical 
analog in soap films, while for parallel orientation nonabelian properties of Wilson loops 
lead to a nonclassical configuration with the same ring but the hole filled by the film. The 
second mechanism arises due to two-gluon exchange between loops and the corresponding 
amplitude is 0(g 4 ). In the confining background and at large N c this simple picture of 
two contours connected by two gluon lines is transformed into a new geometry of two new 
composite loops, as shown in Figs. 4, 5. As a result one has two types of surface configu- 
rations - for small and for large separations between minimal surfaces Si and S2, shown 
respectively in Fig.4 and Fig.5. It is argued in the paper that the configuration generic for 
the scattering corresponds to the Fig.5, and reduces to the (reggeized) glueball exchanges 
between loops while for the case of decay and sea quark loop effects both nonperturba- 
tive and perturbative mechanisms are important with small separation between Si and 
S2. Results obtained in the paper provide a basis for a systematic development both in 
direction of nonperturbative approach to the hadron scattering and in the direction of the 
theory of strong hadron decays. 
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Figure 1: Planar geometry of the Wilson loops for S™ m C S™ m (the case of coinciding 
orientations). 
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Figure 2: Nonperturbative interaction of rectangular Wilson contours. The minimal 
string profile for the case of opposite orientations is given by "enveloping" geometry. 
Leading large area contribution corresponds to annihilation of fluxes along R. 




R 



Figure 3: Nonperturbative interaction of rectangular Wilson contours. The "time-slices" 
of the minimal string profile for the case of parallel orientations are depicted. It is assumed 
that N c = 3. Leading large area contribution corresponds to a single fundamental flux in 
indirection. 
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Figure 4: The same as Figure 5, but with the two-dimensional geometry of Figure 1. 




Figure 5: Dominant two-gluon glueball term responsible for perturbative interaction in 
nonperturbative background at large distances. Gluon propagator lines are replaced by 
double fundamental lines in large N c limit. 
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